arXiv: 1507.04009v2 [astro-ph.EP] 1 May 2016 


Mon. Not. R. Astron. Soc. 000, [l|fT7| (2014) 


Printed 3 May 2016 (MN style file v2.2) 


Chondrule Transport in Protoplanetary Disks 


Aaron Z. Goldberg^’^*, James E. Owen^’^f and Emmanuel Jacquet^’'^ 

^ Canadian Institute for Theoretical Astrophysics, 60 St. George Street, Toronto, MSS 3H8, Canada. 

^Department of Physics and Astronomy, McMaster University, 1280 Main Street W., Hamilton, Ontario, Canada, L8S 4M1. 
^Institute for Advanced Study, Einstein Drive, Princeton NJ, 08540, USA 

‘^Institut de Mineralogie, de Physique des Materiaux et de Cosmochimie, Museum National d’Histoire Naturelle, 

CP52, 57 rue Buffon, 75005 Paris, France. 


3 May 2016 


ABSTRACT 

Chondrule formation remains one of the most elusive early Solar System events. 
Here, we take the novel approach of employing numerical simulations to investigate 
chondrule origin beyond purely cosmochemical methods. We model the transport of 
generically-produced chondrules and dust in a ID viscous protoplanetary disk model, 
in order to constrain the chondrule formation events. For a single formation event we 
are able to match analytical predictions of the memory chondrule and dust popula¬ 
tions retain of each other (complementarity), hnding that a large mass accretion rate 
(> 10“^ Mq yr“^) allows for delays on the order of the disk’s viscous timescale be¬ 
tween chondrule formation and chondrite accretion. Further, we find older disks to be 
severely diminished of chondrules, with accretion rates < 10“® Mq yr“^ for nominal 
parameters. We then characterize the distribution of chondrule origins in both space 
and time, as functions of disk parameters and chondrule formation rates, in runs with 
continuous chondrule formation and both static and evolving disks. Our data suggest 
that these can account for the observed diversity between distinct chondrite classes, 
if some diversity in accretion time is allowed for. 
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1 INTRODUCTION 


The earliest witnesses of the physical and chemical processes 
that took place during the formation of the Solar System, 
4.57 billion years ago, belong to the type of meteorites called 
chondrites. Chondrites are primitive, having mostly escaped 
differentiation processes on their parent asteroids, and are 
composed of various roughly millimeter-sized inclusions na¬ 
tive to the solar protoplanetary disk set in a fine-grained 
matrix. While the oldest among these inclusions are the re¬ 
fractory inclusions, presumably high-temperature conden¬ 
sates from the earliest phases of the disk, the most abundant 
are millimeter-sized silicate spheroids known as chondrules, 
which formed as a result of the solidification of molten 
droplets | Amelin et al.|[20l0 Connolly &: Desdi||2004 1. The 
matrix itself is a composite mixture of micron-size presolar 
components, protoplanetary disk condensates, and possible 
by-products of chondrule-forming events i Huss et al.|2005 |. 

Although theories abound, the nature of the chon¬ 
drule formation events (CFEs) remains very mysterious (e.g. 
|BischofHl998| [Connolly fc Desch|2004[ ) . There are two main 
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categories for CFE models, more or less equally backed in 
the latest literature: “planetary” and “nebular”. The former 
may require phenomena such as collisions between plan- 


etesimals in the disk (Symes et al. 


1998 


fc MovshovitzpOllj [Sanders Scott||2012 [Johnson et al. 


Asphaug, Jutzi 


20151, and a popular example of the latter relies upon shock 


waves passing through the disk (Desch & Connolly 2002 
Morris & Deschp010[ [Morris et b.||2012| |Boley, MoGs fcj 
Desch||2013 1. 


A key constraint is that chondrules, even in single chon¬ 
drites, have ages spread over 0-3 Myr after the condensa¬ 
tion of refractory inclusions. Not only does this suggest that 
CFEs occurred repetitively over the evolution timescale of 
the solar protoplanetary disk, as disks typically last ~ 3 
Myr around young stars ( Williams fc Cieza|2011 |, this also 
indicates that chondrules did not generally accrete immedi¬ 
ately after formation, but spent a few Myr as free-floating 
objects in the gaseous disk. Chondrules may have been sig¬ 


nificantly redistributed by aerodynamic forces (e.g. Cuzzi 
fc Weidenschilling|2006 Jacquet|2014b f At the same time, 


chondrules seem to form distinct populations in different 


chondrite groups, constraining any disk-wide mixing (Jones 
2012 Alexander &: Eb^[2012 b 


Transport of chondrules may also be crucial to un- 
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derstand the compositional diversity of chondrites (Jacquet 


2014b I. Indeed chondrites are not a uniform class, but are 


currently comprised of 14 chemically, isotopically, and pet- 
rographically distinct groups, which to hrst order may be 


subsumed in two superclans (Kallemeyn, Rubin & Wasson 
1996 Warren|2011 1: carbonaceous chondrites (CCs; matrix- 


rich, with mostly solar composition - especially for the 
Cl group - e.g. in terms of the Mg/Si ratio), and non- 
carbonaceous chondrites (including ordinary and enstatite 
chondrites, which are matrix-poor and further from solar). 
One school of thought is that the chondrite spectrum reflects 
varying mixing proportions between a high-temperature 
component dominated by chondrules and a low-temperature 


component believed to be Cl composition matrix (Anders 
1964 Zanda et al.|2006 Zanda, Humayun fc Hewins|2012 1; 
this is the two-component model. Yet this view is challenged 
by evidence that matrix in carbonaceous chondrites is non¬ 
solar, e.g. in terms of its Mg/Si and Fe/Si ratios, and nei¬ 
ther are the coexisting chondrules, although the bulk car¬ 
bonaceous chondrites are close to solar for these elemen¬ 


tal ratios (Bland et al. 2005 Bezel & Palme 20101. This 


matrix-chondrule complementarity - if not due to analyt¬ 


ical artifacts or parent body processes (Zanda, Humayun 


& Hewins 20121 - suggests a genetic relationship between 


the two, with chemical exchange in a closed-system reser¬ 
voir ( Bezel fc Palme||2010 Palme, Bezel fc Ebel|[20l5 L For 
example, if Si evaporated from the chondrules, it could have 
preferentially recondensed on matrix grains owing to their 
larger total surface area. This certainly constrains transport 


of solids since chondrule formation (Jacquet, Gounelle & 
Fromang|2012 1. 


Clearly, the transport of chondrules is important to un¬ 
derstand, not only to better comprehend chondrites as “end 
products” of disk processes, but as a constraint for CFEs 
themselves. Yet, despite considerable amount of work de¬ 
voted to the transport of solid particles in general, as part 
of the effort to understand chondrite accretion (see e.g. re¬ 
views by ICuzzi fc Weidenschilling||2006[ [Chiang fc Youdin] 
2010), the transport of chondrules per se has been little 


investigated specifically in the literature, save as part of 
the - purely analytical - synthesis on chondrite component 


transport by Jacquet, Gounelle & Fromang (2012). In con¬ 


trast, several numerical studies on the transport of refrac- 


tory inclusions already exist (e.g. 

Cuzzi, Davis & Dobrovol- 

skis|2003||Hu|20101|Ciesla|20101|Boss, Alexander & Podolak 

20121 Yang & Ciesla|2012 j. This is 

certainly understandable: 


the origin of chondrules is so ill-understood (even compared 
to refractory inclusions) and the disk properties similarly 
uncertain that such an endeavour may seem premature or 
arbitrary at best. Yet, as time goes by, it does appear that 
chondrule formation and disk transport stand little chance 
to be solved entirely separately, and that a first attempt to 
mutually constrain these interwoven problems in numerical 
simulations has to be made. This is the purpose of this work. 

In this work, we simulate the transport of chondrules, 
chondrule precursors, and dust grains in ID models of “con¬ 
ventional” turbulent gaseous disks and monitor the chon¬ 
drite composition expected at any given time and heliocen¬ 
tric distance to compare with observations. We hrst assume 
a hxed, initial pair of chondrule and dust populations, focus¬ 
ing on the question of matrix-chondrule complementarity. 
We then allow for continuous chondrule formation, follow¬ 


ing simple prescriptions throughout the simulation time, to 
understand the diversity of origins of components (age, he¬ 
liocentric distance of formation) present in individual chon¬ 
drites. 

The outline of this paper is as follows. In Section 
we outline the relevant equations for gas and solid dynam¬ 
ics in the protoplanetary disk, along with an explanation 
of onr numerical methods. Section contains our model of 
investigating the evolution of a single chondrule and dust 
population following a single CFE, and Section extends 
these calculations to disks with multiple CFEs. We discuss 
the implications of our results in Section and, in Section 
we conclude. 


2 DISK MODEL AND METHODS 


We consider an axisymmetric disk in cylindrical coordinates, 
with R the heliocentric distance and 2 the height above the 
midplane. We work within the thin-disk formalism appro¬ 
priate for protoplanetary disks (e.g. Pringle|1981 1 and work 
with vertically averaged quantities. The snrface density of 
gas, dust, or chondrules is defined by 


/ + 00 

pdz, 

-OO 


( 1 ) 


where p[R, z) is the mass density. 


2.1 Gas disk 


At radii R > 0.5 AU, the disk is passively heated (e.g. Chi- 


ang fc Goldreich|1997 D’Alessio, Calvet fc Hartmann|2001 ) 
and as such the temperature (T) proHle is taken to be a time- 
independent power law of heliocentric distance T oc 
with q — 0.5 ( Kenyon fc Hartmann|[r987 ) unless otherwise 
noted. For a vertically isothermal disk the scale height (H) is 
given hy H = Cs/fl, where Cs is the isothermal sound speed 
and n = -yGMe/R ^ is the Keplerian angular velocity. Fol¬ 
lowing Owen (2014) we set the aspect ratio in all of our 


calculations to: 




l/2-q/2 


( 2 ) 


Owing to turbulent angular momentum transport, the 
gas surface density evolves following (|Pringle||1981[) 


asg 

dt 


3 JO 

RM 




( 3 ) 


where the “turbulent viscosity” is given by the a formalism 
(Shakura k, Sunyaev||1973), such that: 


v{R) = aCsH. 


( 4 ) 


For constant a we find o oc R, which is consistent with the 


observational diagnostics of disk evolution (e.g. Hartmann 
et^jdJ1998jJAi^rews et al.|2009p . Following |Oweii, Ercolario 
& Clarke '( [2011 ) we set a = 2.5xl0“®to match the disk life¬ 
times and accretion rates within the X-ray photoevaporation 
model. We note here that, while we have chosen particular 
scales for many of the models we discuss here (steady-disk 
models), our results can be re-scaled to other chosen val¬ 
ues of a and H/R. As such we choose to work with natural 
scales of the disk. 
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Specifically, a natural timescale of the problem is the 
viscous timescale tvis{R) (namely the time over which a gas 
parcel at distance R is accreted by the Sun), given by: 

UUR) = ^ 

= 0.04Myr f ^ 

^ \1A\J J V 0.04 J 
(2.5 X 10-3) ■ 

The strength of the turbulence (and as such the value 
of q) is very uncertain, even at the order-of-magnitude level. 
Therefore, most of the times discussed in the results will be 
normalized to tvis so as to be independent of a (for a given 
accretion rate). 

For timescales longer than tvis{R), the surface density 
should approximate the steady solution with a uniform mass 
accretion rate M = —27ri?EgUij {ur being the gas radial 
velocity), given by: 


E 


g 


M 

3tvu 



( 6 ) 


where i?, is the disk inner edge, and a zero-torque boundary 
condition is applied. Therefore, with our setup at radii S> R* 
the surface density follows E oc R~^. 


2.2 Evolution of solids 


A population of solids, whether chondrules, chondrule pre¬ 
cursors, or dust grains, evolves in the disk due to (i) advec- 
tion by the mean gas flow, (ii) turbulent diffusion, and (iii) 
drift because of finite size leading to partial decoupling with 
the gas. The evolution of the surface density Es of any solid 
population obeys (e.g. Jacquet, Gounelle fc Fromang|2dl2 1: 


dE, 

dt 


i A 

RdR 


R 


"^4 



Es (u_R -f Udrift) 


(7) 

where D = v/Sc is the radial gas diffusion coefficient. Sc 
the radial Schmidt number parametrizing the strength of 
the dust diffusion, henceforth taken to be 1 (unless otherwise 
noted), and 


_TdP 


( 8 ) 


is the solid particle radial drift velocity arising from the gas 
pressure gradient. P{R) = pc^ is the midplane gas pressure 
and 


T = 


TT Patt 
8 pCs 


(9) 


is the stopping time (due to gas drag) of a solid particle, 
with pa and a the particle’s internal density and radius, re¬ 
spectively ( [Jacquet, Gounelle fc Fromang||2012[ ). Through¬ 
out this papeiT^eTake^so'^'OA'gTcni^'^OT (millimeter-size) 
chondrules and aggregates, and psO = 1 x 10“"* g/cm^ for 
(micron-size) dust grains. Note that we work with solid/gas 
ratios ^ 1 so that our results are linear functions of the 
initial solid abundances, and thus relative variations of solid 
abundances are independent of any assumption about their 
initial abundances. 
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We can then dehne the dimensionless gas-solid decou¬ 
pling parameter 


_ rir _ TT PaO 

^ “c7 ^ 2 ^ 


( 10 ) 


as a measure of importance of the drift contribution 
I j Jacquet, Gounelle fc Fromang||2012| ). For a steady disk, 
using Equation (|^, one finds that at R^Rt'. 


S = 


SC Pa aCa 

2 MQ 


0.2 


Pa a 


0.1 g cm“ 


M 


10-8 Mo/yr 


/ (H/R)(lAU) y 
V 0.04 J 


R 

1 AU 


3 / 2-9 


( 11 ) 


which is independent of a for a given value of M. 

No sink term corresponding to chondrite accretion is 
considered here. Nevertheless, at each time and radial lo¬ 
cation, one can define a potential chondrite composition 
that would result from partial accretion of local material. 
We assume that the potential chondrite’s chondrule/matrix 
makeup is representative of the location in question; that 
is, we ignore any possible “accretion bias” as argued by 


Jacquet, Gounelle & Fromang (20121; Jacquet (2014a|. This 


will give us insight into how chondrite composition evolves 
in time and space, and whether there is a spatiotemporal 
“window” matching the observations. 


2.3 Quantifying matrix-chondrule 
complementarity 


As part of the evaluation of the “potential chondrite” com¬ 
position, we will seek to quantify the possible complemen¬ 
tarity between chondrules and matrix as mentioned in the 
introduction. Complementarity cannot simply be a function 
of how close the bulk composition is to solar (that is, close to 
the Cl chondrites, which are deemed to best represent solar 
abundances), for the individual chondrules and matrix may 
also be so close to solar that this is not a real constraint on 
their genetic relationships (not to mention possible analyti¬ 
cal biases or secondary effects; jZanda, Humayun fc Hewinsj 


20121. Thus, we need to factor out the original closeness 


of the different components to solar. Hence, we shall define 
here, for a fiducial chemical element X, a complementarity 
parameter as: 


- _ [A']ch [Ajbulk [A']mx [A^lci 
^ m,,,, - [A]b,ik [X],h-[X]ci’ ^ ' 

where [A'Jch, [A'Jmx, and [Xjbuik are the abundances (by 
mass or normalized to some major element) of element X 
in the chondrules, matrix, and whole-rock of a sample, re¬ 
spectively, with [A] Cl being the Cl value. The motivation 
behind this specific form is that, assuming that there is a 
value rci of the matrix/chondrule ratio r for which a mix 
of the observed chondrules and matrix is solar, reduces to 
rsampie/rci, which is directly measurable in simulations. It 
further does not depend on the chosen initial chondrule/dust 
ratio, nor any chemical assumption in the simplest case of a 
single pair of chondrule and dust populations as investigated 
in Section is of course unity if the bulk chondrite is solar 
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Table 1. Magnesium to silicon ratios for samples of the carbona¬ 
ceous chondrite types CR, CV, CO, and CM, and three non- 
carbonaceous chondrites. All values are normalized to Cl elemen¬ 
tal abundances. Mg/Si is reported for chondrules, matrix, and 
whole-rock values for each chondrite, and used to calculate the 
complementarity parameter (^. Unless the reference mentioned a 
specific Cl Mg/Si value, the Cl value chosen for normalization 


(0.90) was taken from |Palme Sz Jones| ||2005l). ^Hezel & Palme 

1 2 OIOI 1 , ^Klerner et al.|(|2001|| 

"lEbel et al.|(|2008|l, iLodders Sz Fe- 

gley|(|1993t 'IMcSween & Richardson|l 

1977|l, ^Zolensky, Barrett 

& Browningl ||1993ll, IMason 

1196311, ' 

Rubin Sz Wasson] i 

19871, 

^Grossman & Brearleyl (|2005 

, ^^Ahrens, Willis & Erlank|l 

19731, 


^^Huss^trir^OOSll, ^^Grossinan et al.(|T985y7^^R^™b^idi7Hous-| 

|ley &: Rajanj l|1984|l, El Goresy et al.| | |1988[ l, ^ ^Berlin| | |2009[ l, 
^'^Grossman fc Wasson (|1983[l. 


Type 

Name 

Mg/Sici 

(chon¬ 

drule) 

Mg/Sici 

(ma¬ 

trix) 

Mg/Sici 

(whole- 

rock) 

c 

CR 

Renazzo 

1.15^’^’^ 

0.71'^’^’^’ 

’ 1.01'' 

0.93 

CV 

Vigarano 

1.18® 


1.03'^ 

0.70 

CO 

Kainsaz 

I.I 91 

0.85i'®-“ 

0 . 991 “ 

1.07 

CM 

El-Quss 

Abu 

Said 

I. 23 I 

O. 73 I 

l.OOi’ii 

0.99 

E 

Qingzhen 

0 . 941 ^ 

O. 751 ® 

0.831'! 

-5.90 

K 

Kakangari 

1.02^® 

I.O 41 ® 

0.981® 

1.52 

0 

various 

O. 93 I® 

10 

00 

l> 

0 

0.88i“ 

-1.69 


in X. As such, the closer ^ is to 1, the more complementary 
a sample is. 

In order to get a feel for real values of (^, we take the 
Mg/Si ratio as our [X] (but the calculations herein are not 
unique to this specific chemical parameter and thus should 
apply to any claimed complementary relationship). Com¬ 
positional data for a number of chondrite classes are tabu¬ 
lated in Table |2.3| We see that CCs tend to have ^ within 
~ 10% of 1, with the CV chondrites deviating most from 
this G ~ 0.7). In contrast, enstatite and ordinary chondrites 
have negative values because both their chondrules and ma¬ 
trices have a subsolar Mg/Si ratio (so that no combination 
thereof can restore a Cl chondritic composition); again com¬ 
plementarity is more a property of carbonaceous chondrites 
than of chondrites in general. In this paper, unless otherwise 
noted, we will thus adopt the convention that a chondrite is 
“complementary” if its Q lies between 0.7 and 1.3. 


2.4 Numerical methods 


In order to solve for the transport of solids and any possi¬ 
ble chemical mixing in a time-dependent manner we must 
proceed numerically. Our numerical method integrates the 
evolution of the gas Q and dust Q surface densi¬ 
ties explicitly on a staggered non-uniform mesh using fi¬ 
nite volume operators. It is identical to that described in 
[o^poTit and has previously been used for astrophysi- 
cal and cosmochemical applications by [Owen fc Armitage| 
(20141; Owen & Jacquet ( 2015[ ). The method is second- 
order-accurate in space and first-order-accurate in time. For 
the advection term in Equation ([^ we use second-order re¬ 
constructions and Van-Leer limiters at cell boundaries. At 
the inner boundary we adopt zero torque boundary con¬ 
ditions for the gas and free-outflow boundary conditions 


for the solids. In the simulations with a steady disk pro¬ 
file we have constant mass-flux boundary conditions at the 
outer edge, in both the gas and the dust components. In 
the simulations with an evolving disk we apply at the outer 
edge a zero-torque boundary for the gas and free-outflow for 
the solids, but we note that in these simulations the outer 
boundary is chosen to be at a distance such that it does not 
affect the evolution. 

Unless otherwise noted, the simulations were run with 
800 grid elements spaced uniformly in 7?^^® between i?in = 
3 X 10^^ cm and Rout = 3 x 10^^ cm. 


3 EVOLUTION OF A SINGLE PAIR OF 

CHONDRULE AND DUST POPULATIONS 

In this section, we investigate steady disk models with a 
fixed, initial pair of chondrule and dust populations, that is 
only consider a single CFE. 


3.1 Complementary start 


We first consider the case where chondrules and dust are 
cogenetic and originally in complementary proportions. We 
assume that the initial populations are within a “top hat” 
of width L centred at radius Rcentre- Within the limits 
R = Rcentre “ L/2 and R — Rcentre + 7//2, the chondrule 
and dust surface densities have Ech = Ed = 0.01 • Eg (but 
note that the complementarity parameter C, at any time is 
independent of the initial solid/gas ratios chosen). We evolve 
the populations over 2Myr for values of Rcentre, L/Rcentre, 
and M ranging between 0.3-30 AU, 0.05-1.5, and 10“^°- 
10“® M0yr“^, respectively. 

Over time, the chondrule population tends to drift in¬ 
ward and its width grows larger than L, exemplified by Fig¬ 
ure [Tlr. Its inward drift being faster than that of the dust, 
decoupling between the two originally complementary pop¬ 
ulations progressively occurs. To make this quantitative, we 
calculated the fraction w (the “complementary fraction”) of 
the chondrule mass within the simulation domain for which 
the complementarity parameter (here independent of any 
assumption on the composition of chondrules relative to ma¬ 
trix) is in the conventional complementary range (0.7-1.3; 


see Section 2.31. By construction, the system is completely 


complementary {w = 1) at time t = 0, but w monotoni- 
cally decreases with time (Figure]^). We call thaif the time 
at which w = 0.5, which serves as a measure of the chon- 
drule/dust decoherence timescale ( [Jacquet, Gounelle fc Fro 
man 


In Figurej^ we plot thaif/tvis(Rcentre) versus SIR centre) 
(evaluated for chondrules). Under this normalization, the 
points for a given L/Rcentre value collapse on a single line. 
This is due to the fact that the solid evolution equations 
([^ can be non-dimensionalized in terms of only R/Rcentre, 
t/tvis (Rccntre) , and S{R centre ) for a steady disk. Figures 
and show how this behaviour is modified if the threshold 
complementary fraction or the range of g deemed comple¬ 
mentary are changed, respectively. 

Several regimes can be distinguished depending on the 
value of S'. 

For very large S, specifically S > Rcentre/R, gas-solid 
drift is dominant and diffusion hardly affects the chondrule 
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Figure 1. a. Non-dimensionalized time i/tvis(-Rcentre) versus complementary fraction w, for .Rcentre = 0.75AU, L/Rcentve = 0.5, and 
M = 1.36 X 10“® M 0 yr“^. The mass fraction of complementary chondrules is seen to decrease monotonically from 1, with fhaif ~ 2fvis- 
b. Space-time plot of the chondrule population evolution for the same parameters. The plot is coloured according to chondrule mass 
fraction, overlain on top of which are contours of constant C,. Both plots have a common y-axis, with non-dimensionalized time increasing 
from top to bottom on the graphs. 




Figure 2. Log-log plot of t{w = /)/ivis(7?centre) versus S for 
varying -Rcentre and M, evaluated at .Rcentre, M^ and various val¬ 
ues of /. The value of w drops below larger values of / in shorter 
amounts of time, but the slopes do not rely heavily on the chosen 
value of / for any value of S. 


Figure 3. Plot of t^aif various definitions of complementarity, 
all with L/Rcentre = 0.1. Values of fhaif are normalized to the 
value corresponding to the standard complementary range (0.7 ^ 
C, ^ 1.3) in the paper. 


and dust profiles, which are essentially translated inward, 
t/; = 0.5 then corresponds to the case where the chondrule 
population is shifted by I//2 relative to the dust population, 
which occurs for the time 


thalf — 


L/2 

I't’driftI 


Lp dP 
2r dR 


L 

2 


tvis (R-ci 


S(Rc, 


ainP 

ainR 


(13) 


(neglecting the drift of the dust grains relative to the gas), 


yielding (for a steady, constant a disk; i.e., P oc ^1"^^ 

, /p^ = (ir^) (6-g)-'5(R.onne)-V (14) 

^vis t-^centre J \-^centre / 


For 1 ^ S' i?centre/r, diffusioual widening of the ini¬ 
tial populations cannot be ignored, and in fact rapidly dom¬ 
inates over the initial width L, delaying the decoherence be¬ 


tween chondrules and dust (Jacquet, Gounelle & Fromang 
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2012). We show in appendix [A| that 


^half — 


2D 


2.2972 

This leads to the relation 


■ dP/dR 


(15) 


^half 


^vis {R 

Cl 


2.2972 ScS(R.. 
0.385'c 


■centre 

1 


/91nPy 

)2 


((3-g/2) 5(R centre)) 


(16) 


For S <C 1, both chondrules and dust should follow the 
gas and be advected to the disk’s centre by time ~ tvis- How¬ 
ever, complementarity may linger for even longer timescales 
while part of the chondrules and dust are diffused outward 
(as w is indeed defined in terms of the remaining chondrule 
population). We may estimate that complementarity will be 
lost when diffusion will have reached the Sr = S ■ Sc = 1 
line, which is the cutoff for outward diffusion for the chon¬ 


drules but not for the dust (Jacquet, Gounelle & Fromang 


2012), that is on a timescale 


thalf ^ RisISr — 1) — tvis (Rcentre 


\ -^centre J 


(17) 


and so 

^half 

^vis (7?cen1 


= (*S'(-Rcentre ) • 5 c)1-2/(3/2-9) . 


(18) 

For our q value, this corresponds to a S~^ dependence as 
observed. Additional runs with alternative q values change 
the slope of the log-log plot in Figure in accordance with 
this prediction. Clearly, for 5^1, complementarity can be 
maintained for timescales as large as the viscous timescale, 
as advocated by [Jacquet, Gounelle fc Fromang| ( |2012[ ). 

We have re-run the trials with only two M values, a 
wider range of i?centre values, and 16 times the spatial reso¬ 
lution, to validate convergence to the analytical results (Fig¬ 
ure]^. The —1 slope from (18) and exact solution taking 
q = 0.5 for (16) robustly match the data, but the predicted 
behaviour 01 at large S is oft by a factor of about 2/3. 
The latter may be attributed to rapid diffusion at the sharp 
top hat edges with which the populations are initialised, as 
it improved slowly with increased spatial resolution. Setting 
Sc = 10® to nullify the effects of diffusion for a few trials 
with 5 1 resulted in agreement with the analytical pre¬ 

dictions to within 2%. 


3.2 “Two-component” start 

As a comparison, we have run simulations with spatially sep¬ 
arated chondrules and dust populations, with the former and 
the latter initialized within the inner and outer halves of the 
top hat, respectively, with L/Rcentre = 1 and 1.5. This mim¬ 
ics a “two-component” picture where chondrules are derived 
from the inner disk, whereas the matrix represents pristine 
outer disk matter. Here, chondrules and dust are assumed to 
initially be in “global” complementary proportion^ Since 

^ It should be noted that this may not exactly reflect the two- 
component model discussed in cosmochemistry, since the latter 


they are spatially distinct, the initial complementary frac¬ 
tion w as defined above is zero. With the passage of time, 
chondrule particles diffuse outward and some dust particles 
drift inward, thus w increases and reaches a maximum of 
w = 0.5-0.6 at time t ^ thaif, then decreases back to 0 as 
chondrules drift past the inner boundary (e.g. Figure [^. 
Regardless of the values of M and Rcentre, the total amount 
of time with w ^ 0.5 never exceeds 2.3 x 10“®tvis(Rcentre). 
Complementarity, provided it is indeed a pristine feature of 
carbonaceous chondrites, is therefore much more difficult to 
obtain in a two-component picture vis-a-vis a single reservoir 
model. 


4 CONTINUOUS CHONDRULE FORMATION 

While the previous section considered a single initial CFE 
and the evolution of its products, we now investigate mod¬ 
els with multiple epochs of chondrule formation. Following 
the results of the previous section, we shall consider that 
chondrule and matrix grain formation were not spatially 
distinguished, save for whichever fraction of an initial Cl 
chondritic matrix component escaped CFEs until chondrite 
agglomeration. We first present our prescriptions for chon¬ 
drule formation, then simulations for steady disks, and fi¬ 
nally results for an evolving disk scenario. 


4.1 Initial conditions and chondrule formation 
model 

We consider initially two populations of Cl chondritic com¬ 
position (in terms of a fiducial element X) particles, one 
micron-sized (the dust) and another millimeter-sized (“ag¬ 
gregates”). Neither growth nor fragmentation are modelled. 
The solids’ initial surface density profiles follow the static 
(for a static gas disk) derived by 


%^(-R) = exp(/(R)) r dR' (19) 

with € = 0.01 and 

I{R) = ^iR')dR', (20) 

with the outer boundary conditions for both popu¬ 
lations set to maintain the accretion rate (neglecting 
the diffusive flux) 

Afd.agg = 27r77out Ed,agg (Rout )u_R;d,agg (Rout) “ eA/(Rout). 

( 21 ) 

Chondrule formation is modelled with a chondrule pro¬ 
duction function (g(R)) defined such that, for any time in¬ 
terval dt, a fraction gdt of the millimetre-sized bodies (orig¬ 
inally only aggregates, but later also including previously 
formed chondrules) at any given location is converted into 


solution Equation A.6 


Jacquet ^ Robert (2013 


often assumes the matrix to be solar (e.g. |Zanda, Humayun &:| 
|Hewins||2012| |, although |GonzaIez| ( |20l'^ suggests Cl chondrites 
(which are used to represent solar abundances) are systematically 
volatile-enriched relative to the Sun. Here, our point is merely to 
evaluate how difficult complementarity is to obtain from non- 
cogenetic matrix and chondrules. 
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Figure 4. Log-log plot of thalf/^via versus S for varying -Rcentre, M, and L/Rcentre- Both are evaluated at -Rcentre and M, and overplotted 
in solid blue and red lines are the analytical predictions. Red and blue points correspond to L/Rcentre = 0.05 and L/Rcentre = 0.1, 
respectively; plus signs and circles correspond to M = 1.85 X lO“^®M 0 yr“'^ and M = 2.15 X lO“®M 0 yr“^, respectively. The blue and 
red points appear as one red point for S 1 , as they completely overlap. 


new chondrule^ The same fraction of the existing dust be¬ 
comes a new dust population. As a proxy for chemical ex¬ 
change between dust and chondrules, the concentration of X 
in the new chondrules and dust is fixed such that the local 
chondrule/dust “partition coefficient” 

( 22 ) 

is a constant taken to be 1.37, an average calculated for 
CC data, with the additional constraint of conservation of 
bulk chemical abundance during the CFE. When assessing 
potential chondrite composition, the remaining primordial 
millimeter-sized aggregates at the location and time of con¬ 
sideration were counted among the matrix component. 

We chose g to be zero outside of a prescribed chondrule 
forming region (CFR) between 0.5 AU and Rcfe (which was 
varied from 3 to 25 AU) and to be of the form 

" tvis(lAU) ^ (l^) 


^ Other prescriptions were tried where part of the dust is also 
converted into chondrules but were not found to significantly dif¬ 
fer from our standard “conservative” prescription. For the sake of 
simplicity, given the arbitrariness of the chondrule recipe anyway, 
we shall solely focus on the latter. 


in the CFR. A is a constant, corresponding to the fraction 
converted at 1 AU within one viscous timescale, and S indi¬ 
cates the radial dependence of chondrule formation, ranging 
from 0 (no dependence) to 2 . 

Numerically, the CFR was divided into several radial 
bins, and the CFEs were discretized to take place with a 
finite period dt. For each CFE time and each radial bin, a 
(Eulerian) chondrule and a dust population (originally con¬ 
fined to the radial bin in question) are followed, that is are 
first created at the CEE time in question (with subtraction 
of the corresponding precursors from the previous popula¬ 
tions), and then have their radial distributions updated for 
the remainder of the simulation following Equation 0 - For 
simplicity, since information from the exact provenance in¬ 
side a radial bin is lost, a single (average) composition was 
assigned to each chondrule and dust population, correspond¬ 
ing to the average one resulting from Equation (22 1 and mass 
conservation. 


4.2 Steady disk simulations 

We have run simulations using steady gas disks with M from 
lO~®M 0 yr“^ to lO~^M 0 yr“^ for 0.84Myr ~ 21tvis(lAU). 
We varied S (0 and 2), Rcfe (3, 10, and 25 AU), and A 
(0.075 and 0.75). 12 radial bins were used with the CFE 
period set at dt = O.OlMyr. 
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Radius (AU) 


Figure 5. a. Non-dimensionalized time t/tvis(^centre) versus complementary fraction w, for -Rcentre = 4.1 AU, L/Rcentie = 1-5, and 
M = 1 X lO“^M 0 yr“^. The mass fraction of complementary chondrules is seen to steadily increase from in = 0 to tu Ri 0.25 before 
decreasing back to 0, with no t{w ^ 0.5). b. Space-time plot of the chondrule population evolution for the same parameters. This plot 
is coloured according to chondrule mass fraction, overlain on top of which are contours of constant i^. Both plots have a common y-axis, 
with non-dimensionalized time increasing from top to bottom of the graphs. 



Figure 6. Space-time plot of the chondrule-to-dust surface density ratio, overplotted with contours of constant a. Cl-normalized bulk 
composition and b. complementarity parameter Q. This simulation has M = 10“® MQyr“^, 5 = 0, i?CFB = 25 AU, and A = 0.75. 


Figures and show space-time diagrams of the 

chondrule/dust ratio, with contours of the bulk composi¬ 
tion as well as the complementarity parameter, for the case 
A = 0.75 (0.075 in Figure]^, <5 = 0, Rcfe = 25 AU, and 
M = lO“®M 0 yr“^ (lO“^M 0 yr“^ in Figure]^. We wit¬ 
ness the initial increase of the chondrule abundance (initially 
zero) until some steady-state is reached with their loss by 


advection and drift into the Sun. A radial gradient of the 
chondrule/dust ratio sets in, in part due to the limited ex¬ 
tent of the CFR, but here mainly controlled by the differen¬ 
tial drift between chondrules and dust (more or less coupled 
with the gas). As we mentioned previously, this starts to 
have an effect when S' ~ 1 for the chondrules, corresponding 
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Figure 7. Same as Figure]^ but with a lower A = 0.075. Sch/^d decreases accordingly by one order of magnitude, Cl-normalized bulk 
composition remains within 10 % of 1 at all times, and the complementary range of ^ moves to smaller radii (and remains constant). 



Figure 8 . Same as Figure]^ but with a higher M = 10 ^M 0 yr In contrast to Figuresand [ t] chondrules occupy a portion of the 
disk out to 20 AU. 


to (using Equation (111): 


^ VlO-8M0/yri?/i?(lAU) 

(24) 

Indeed, as M decreases (compare e.g. Figures to [^, the 
chondrules become more and more concentrated in the in¬ 
ner disk. Compositionally, this results in a more chondrule- 
enriched innermost region. Beyond this in the CFR is a 
chondrule-depleted region, which shows a slightly subsolar 


/ 2 \ 1 /{ 3 / 2 - 9 ) 

O.lg/cm \ 

Psa J 


composition in term of our fiducial “chondrule-loving” ele¬ 
ment X for the lowest M = 10~® MQ/yr, with yet further 
out a chondritic composition dominated by the unprocessed 
dust/aggregates. Only for the highest M = 10“^ M 0 /yr 
(Figure]^ do we see complementarity linger for a few fvis, 
in accordance with the previous section. 

For smaller A values, the steady-state chondrule/dust 
ratio is lower and the departures from solar compositions 
are thus smaller. A very rough quantification of the pro- 
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Origin Radius (AU) Time Since Formation /t. (1 AU) 


Figure 9. Origin distributions of chondrules found at i? = 3 AU, in terms of mass fraction, after a. & b. 0.035 Myr (0.9tvis(l AU)), 
c. &; d. 0.405 Myr (10.2tvis(l AU)), and e. &; f. 0.795 Myr (20.0tvis(l AU)), each of which are halfway between CFEs. Here, M = 
10“® M 0 yr“^, (5 = 0, Rcfe = 25 AU, and A = 0.75 are used. Note that the age is defined relative to accretion time, not to the start 
of the simulation. Origin radius is taken as the middle grid cell of each chondrule formation bin. Chondrules are predominantly seen to 
have formed at radii close to their current radius, with an emphasis towards larger radii, and their age distribution falls off exponentially 
with time since formation. These distributions achieve a steady state by £ 5^ 10£vis(l AU). 


duction/loss balance is to consider that the net production 
of chondrules 5 'Sdust+agg in the inner disk is balanced by 
the loss —Ech(l + 5')/£vis (advection+drift) in steady state, 
hence Ech/Sdust+agg ~ gUis/A + S) = /{1 + S) 

(compare Figures]^ and [^. 

Figure shows the age and location of origin of chon¬ 
drules found at 3AU after 0.04, 0.4, and 0.8 Myr (midway 
between the two closest CFEs). From these histograms, we 
see the intnitive result that most of the chondrules at a given 
radius were formed in spatial bins near that radius, with a 
greater emphasis on the bins beyond it due to inward advec- 
tion. The ages of the chondrules at this radius are mostly 
within tvis(3AU) « 0.12 Myr, with a trend toward recent 
formation. The distributions are nearly identical at the later 
two times. This implies that the continuous CFEs lead to a 
constant distribution of spatiotemporal origin of chondrules 
at a given radius within 0.4Myr « 3tvis(3AU), due to the 
continuous replenishment of chondrule material. Let us now 
study how the spatial and temporal sources of these depend 
on model parameters. 

Figure [TO^ shows a plot of the average chondrule age as 
a function of R (i.e., computed at each radial point) mid¬ 
way between the two final CFEs, for M — lO~®M 0 yr“^. 
The average age is seen to be minimum within the CFR, 
which reflects the dominance of locally, and generally re¬ 
cently produced, chondrules. It steadily increases beyond 
Rcfe, owing to the increasing transport time needed, until 
the limit of outward diffusion, where the minute amount of 
particles are homogenized in terms of spatiotemporal ori¬ 
gin due to numerical effects. The 5 = 2 case shows a more 
complicated behaviour, with two local minima, presumably 
because chondrules are essentially produced at small radii. 


amounting to a small effective outer boundary of the CFE. 
It shows an anticipated increase in age relative to the 5 = 0 
case, yet the few chondrules produced at large distances 
dominate further out and cause another “local production” 
minimum. 


Figure [lOt r shows the average source radius of chon¬ 
drules as a function of radial location. Inside the inner 
boundary of the CFR, it shows a plateau value correspond¬ 
ing to chondrules advected dominantly from the inner re¬ 
gions of the CFR, before rising in the CFR, being com¬ 
parable to, but larger than, the local radius, and then 
plateaus again outside of the CFR at R ~ 0.85Rcfe. For 
Rcfe = 25 AU, <5 = 2 shows smaller average source radii 
than 5 = 0 , as expected from the greater concentration of 
chondrule production at shorter heliocentric distances. The 
plateau outside the CFR may be understood as follows: in 
steady state, the distribution of the integrated population of 
chondrules produced at a given heliocentric distance Rsource, 
which has no source (or sink) term there and zero value 


at infinity, is given by Ei/Eg(R) oc ^exp(—llS'/4)/R^/^) 
( Jacquet, Gounelle fc Fromang||2012 1 . Since the right-hand 
side is a fixed function of R, independent of Rsource, the rela¬ 
tive proportions of two integrated chondrule populations do 
not depend on R and merely reflect the relative production 
rates as well as the relative ease of outward transport at the 
sources. The standard deviation of the source radius (Figure 


10:) decreases from about Rcfe/ 4 inwards of the CFR to 


< 1 AU outwards, reflecting that chondrules produced close 
to the Sun are not transported efficiently outside the CFR 
(e.g., because the CFR ends beyond the S = 1 line). 


Decreasing A from 0.75 to 0.075 increases somewhat 
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Figure 10. Chondrule origins after 0.835 Myr, plotted versus R, with M restricted to 10“® M 0 yr“^. a. Average age normalized to tvia is 
plotted. Age is seen to decrease until the end of the CFR, then rapidly increase until the limit of outward diffusion. For 5 = 2, chondrule 
age grows through the middle and end of the CFR. Increasing A simply reduces average chondrule age. b. Average source radius (AU) is 
plotted. This is constant for small R, increases with R across the CFR, then plateaus. Increasing Rcfe increases average source radius, 
and decreasing A increases source radius by a constant factor for radii smaller than the CFR. c. Standard deviation of source radius 
(AU) is plotted. The width of the radial distribution is seen to be constant for small R, decrease across the CFR, and then plateau, for 
(5 = 0. For (5 = 2, the distribution of source radii is more narrow at small R. 


the average ages and radii, as expected from lower degrees 
of recycling (allowing longer average travels), but the effects 
are rather limited (order unity or less) especially at large 
heliocentric distances. This is because with such A values ^ 
1 we are mostly in the g <C tT regime where the behaviour 
of chondrules is dominated by transport, so that the age is 
rather limited by the viscous timescale of accretion to the 
Sun, and variations in A do not affect the relative production 
efficiencies of the different radial locations. Of course, it may 
be noted, in the opposite regime of efficient recycling g ^ 
^visi chondrules, which would be dominantly local, would 
have an average age asymptoting to g~^ oc A~^, with an 

exponential distribution. _ 

Decreasing M (Figure 111 tends to reduce the ages 
(down to half the numerical period for M = 10“® M 0 /yr). 
This reflects weaker coupling of the chondrules with the gas 
{S oc M~^), incurring more rapid radial drift to the Sun. 
The radial standard deviations decrease somewhat, as pre¬ 
dicted by Jacquet, Gounelle Fromang (20121, though this 
is a limited effect. 


4.3 Evolving disk models 

In our final set of calculations, we consider runs where the 
gas disk evolves. The surface density profile was initialized 
so as to match the zero time Lynden-Bell & Pringle (19741 
similarity solution 


S(i?,) 


Afdisk (0) 

2'kRRi 


exp 


_R 

Ri 


(25) 


fbfdisk (0) = 0.07 M 0 and Ri = 18 AU, corresponding to the 
disk’s initial mass and scale radius respectively, were chosen 
so as to best fit observed accretion rates and disk fractions 


to the Lynden-Bell Sz Pringle solution (Owen, Ercolano &z 


Clarke 20111. These trials were run for 3.5 Myr (that is, 


the time at which the disk begins to clear; |Owen, Ercolano] 


& Clarke 20111. They had the same range of 5 and i?cFB 


values, with 10 radial bins but with a larger dt = 0.175 Myr, 
since the number of chondrule/dust populations that can be 
handled by the code is limited by the available memory. The 
A values explored thus had to be lower (0.01-0.04) than in 
the steady-state runs to avoid over-processing chondrules in 
single CFEs. 


Figure shows an example space-time diagram of the 
results of this simulation, using the same values of M, Rcfe, 
and S as in Figureand the largest value of A. Even though 
CFEs constantly occur, the area containing the majority 
of chondrule material continuously moves to smaller radii 
(along with contours of constant bulk composition and con¬ 
stant C) as time progresses. This corresponds to the decreas¬ 
ing mass accretion rate onto the central star, explaining why 
exceedingly few chondrules are present at later times (after 
~ 2 Myr, i.e. for M < 5 x 10“® M 0 /yr), in accordance with 
observations from the steady disk trials. 


Figure displays the histograms of source radii and 
epochs at 3 AU at four different times. The histograms are 
fairly comparable, at first glance, to their steady-state coun¬ 
terparts (see in particular the t=1.14 Myr panel, correspond¬ 
ing to M — 10“® M 0 /yr, to be compared with Figure]^. 
There is relatively little evolution in the first 2 Myr, in the 
direction of narrowing in agreement with trends for decreas- 
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Figure 11. Chondrule origins after 0.835 Myr, plotted versus i?, with A restricted to 0.75. a. Average age normalized to tvis is plotted. 
Increasing M shifts the behaviour of average age versus R to higher average age. For low M, average age achieves the minimum possible 
value and plateaus for sufficiently large R. b. Average source radius (AU) is plotted. Decreasing M similarly shifts the behaviour of 
average source radius versus R to shorter values, c. Standard deviation of source radius (AU) is plotted. M does not seem to greatly 
affect the width of the source radius distributions at small R, and significantly decreases it as R approaches ~ 30 AU. 



Figure 12. Space-time plot of the chondrule-to-dust surface density ratio, overplotted with contours of a. constant Cl-normalized bulk 
composition and b. constant C,. This simulation had 5 = 0, .Rcfe = 25 AU, A = 0.04, and the conservative prescription. Chondruies 
occupy a large portion of the inner disk, and the outer radius of the area they occupy decreases as time progresses. Contours with 
near-solar bulk composition and C, within the complementary range robustly overlap with the high Sch/^d region. After ~ 2 Myr, M is 
sufficiently small such that chondruies do not remain in the disk for very long. 
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ing M discussed in the previous subsection. This can be 
seen more generally in Figures and which display the 
average ages and source radii as functions of heliocentric 
distance for the different simulations (similar to Figures 10 
and in the steady-state case) at the same four times as 
Figure 1 13) We note a downturn in average source radius fur¬ 
ther outward, which may reflect a significant contribution 
of chondrules produced close to the Sun transported there 
by the original expansion of the disk, an effect seen by|Yang| 


&: Ciesla (20121 in the context of refractory inclusions. This 


may explain the somewhat shorter source radii and older 
ages overall relative to the steady-state runs. This contribu¬ 
tion of early-transported chondrules becomes dominant at 
3 Myr when newly produced chondrules are rapidly lost, as 
may be seen in the final panel of Figure |14| although all 
chondrules in general are scarce at that point. There is also 
no dip in age in Figure [^unlike the steady-state runs. This 
is presumably partly an artifact of a larger dt, but may also 
relate to the inflow of the early-transported chondrules from 
further out. 


5 COSMOCHEMICAL IMPLICATIONS 

5.1 Matrix-chondrule relationships 

One first conclusion of the above study (in particular Section 
is that a complementary relationship between oogenetic 
chondrules and dust can be preserved for long timescales 
provided that the decoupling between chondrules and gas 
is limited, that is S 1. This should hold in the inner 
solar system for relatively high mass accretion rates (M > 
10 “^ MQ/yr), viz., early in the disk evolution as our contin¬ 
uous chondrule formation simulations (Section]^ verify. We 
also note that low mass accretion rates M < 10“® MQ/yr 
would likely not allow chondrules to be retained efficiently in 
the disk because of radial drift, depending on their produc¬ 
tion rate. Thus evidence of matrix-chondrule complemen¬ 
tarity in carbonaceous chondrites (e.g. |Hezel fc Palme|20ld{ 
Palme, Hezel &: Ebel|[2015( is not inconsistent with signifi¬ 


cant transport between chondrule formation and chondrite 
accretion ( [Jacquet, Gounelle Fromang|2012p . There is fur¬ 
ther no inconsistency with evidence for the presence of a Cl 


chondritic component in carbonaceous chondrites (Anders 


|1964[ |Zanda et al.|20d^ |Zanda, Humayun fc Hewins|2012f 7 

as exemplified by our modelling of primordial Cl dust that 
has been mixed with processed matrix components. This 
does not however prejudge the possibility that the empir¬ 
ical evidence of matrix-chondrule complementarity found 
in carbonaceous chondrites, e.g. as to the Mg/Si ratio, 
is compromised by parent body alteration or instrumental 
biases ( Zanda, Humayun fc Hewins||2012 (. 

Given that higher levels of gas-solid decoupling (higher 
S) do lead to chondrule/dust fractionation, it is then tempt¬ 
ing to associate the nonsolar compositions of the noncar- 
bonaceous chondrites with this effect, which was one of the 


suggestions by Jacquet, Gounelle & Fromang (20121. This 
would, however, be difficult for the Mg/Si ratio. Indeed, 
non-carbonaceous chondrites display a subsolar Mg/Si ra¬ 
tio, while the higher-Mg/Si chondrite components are the 
chondrules themselves. Our simulations suggest that the in¬ 
nermost regions of the disk would be enriched in chondrules 


and give rise to chondrites with an enhanced Mg/Si ra¬ 
tio, unlike observations. True, a correspondingly low-Mg/Si 
region may appear further out, but the effect would be com¬ 
paratively limited, and would predict that non-carbonaceous 
chondrites are depleted in chondrules relative to carbona¬ 
ceous chondrites, the contrary of observations. The trans¬ 
port of chondrules by itself is thus unlikely to have generated 
the difference between carbonaceous and non-carbonaceous 
chondrites, although it could explain part of the diversity in¬ 
ternal to non-carbonaceous chondrites, for example regard¬ 
ing metal/silicate fractionation ( Zanda et al.|20d6 Vernazza 
et al. 20141. A remaining possibility, independent of chon 


drule formation, to be investigated in the future, is the loss 
of amoeboid olivine aggregates (a class of refractory inclu¬ 
sions; Larimer fc Wasson||1988 Jacquet|2014b (. 

We have ignored settling effects in case accretion takes 
place preferentially at the midplane. While this would incur 
little change for S 1 (and thus not affect the comple¬ 
mentarity argument above), this would increase the chon¬ 
drule/dust ratio in chondrites above the surface density ra¬ 


tio (Jacquet, Gounelle & Fromang 20121. But, if accretion 


is significant, the disk, becoming preferentially depleted in 
chondrules, would become richer in dust, which might invert 
the trend. Another possibility is that dust accreted on chon¬ 


drules as rims before final agglomeration (Metzler, Bischoff 


fc Stoeffler|1992 Ormel, Guzzi fc Tielens|2008 1 , so that rep- 

resentativity of chondrites vis-a-vis the whole nebular reser¬ 
voir would be more faithful. 


5.2 Space-time distribution of chondrules 

Taken at face value, our simulations (with a = 0.0025) indi¬ 
cate that the source regions of chondrules of a given chon¬ 
drite would be dominantly in the neighbourhood of the 
chondrite’s accretion region, with a radius standard de¬ 
viation equal to a significant fraction of the overall width 
of the GFR in the inner disk and smaller further out. One 
could accommodate a few distinct contemporaneous chon¬ 
drite groups along the GFR, but certainly not > 14 as ob¬ 
served ( [Jones poT^ . But since the typical age of chondrules 
would be of order tvis (here <C 1 Myr), there would be ample 
room for some of the chondrite diversity to reflect diversity 
in chondrite accretion time (in addition to heliocentric dis¬ 
tance). We further note that the narrow age range would 
be dictated by disk dynamics and, if confirmed empirically, 
would thus not imply immediate chondrite accretion follow¬ 


ing chondrule formation (contra Alexander & Ebel (2012|). 

Yet the radiochronological data suggest a wider age 
range (a few Myr) for chondrules in individual chondrites 
( Villeneuve, Ghaussidon fc Libourel] 2009 Gonnelly et al. 
20121 , although the possibility of secondary disturbances 
is not yet ruled out ( [Alexander fc Ebel||2012] ). This might 
simply indicate that the a value used (0.0025) is too high 
and that the viscous timescale is more comparable to a few 
{oi < 10“^), as might be expected in a dead zone cur¬ 
rently believed to encompass a large fraction of the planet¬ 
forming region of protoplanetary disks ( Gammie|1996 1 . (Re¬ 
call that changing the a for a given mass accretion rate 
merely amounts to rescaling the time so our results ex¬ 
pressed in terms of tvia remain valid.) It could well be that 
a steady state in terms of chondrule relative radius-age dis¬ 
tribution was not yet reached, so that more spatial diver- 
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Figure 13. Distributions of chondrule origins at R = 3 AU, in terms of mass fraction, after a. & b. £ = 0.437Myr (11.0£vis(l AU)), 
c. & d. £ = 1.14 Myr (28.6£vis(l AU)), e. &; f. £ = 2.01 Myr (50.6£vis(l A-U)), and g. &; h. £ = 3.06 Myr (77.0£vis(l AU)), each of which are 
0.05 Myr after a CFE. Origin radius is taken as the middle grid cell of each chondrule formation bin. Chondrules are seen to have mostly 
been formed near their current radius, with a greater emphasis towards smaller radii at larger times. There are chondrules present at 
3 AU that have been formed at all times after 1.5 Myr, with an emphasis towards recent formation. 
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Figure 14. Average chondrule age normalized to £vis(l AU) plotted versus R. Each plot occurs midway between the two closest CFEs, 
with a. t = 0.437Myr (11.0£vis(l AU)), b. £ = 1.14Myr (28.6£vis(l AU)), c. t = 2.01Myr (50.6£vis(l AU)), and d. £ = 3.06Myr 
(77.0£vis(l AU)). Average age is smallest and constant at small R, before smoothly increasing; the onset of this increase correlates with 
the age of the disk. The value of A is seen to be negligibly affect age at most times. In a. and b., an increasing value of Rcfe leads to 
a more gradual increase in average age, and the 5 = 2 case coincides with the Rcfe = 3AU distribution. In d., it is observed that few 
young chondrules are seen at radii beyond 30-50 AU. 
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Figure 15. Average chondrule source radius plotted versus R at a. t = 0.437Myr (11.0tvis(l AU)), h. t = 1.14Myr (28.6tvis(l AU)), c. 
t = 2.01 Myr (50.6tvis(l AU)), and d. t = 3.06 Myr (77.0tvis(l AU)). For all times, average source radius is approximately the same and 
constant over small radii. For the first three times, average source radius increases smoothly across the CFR, with the increase beginning 
at a larger radius for later times. In d., source radius is constant up to i? = 30-40 AU, where it sharply increases to some average CFR 
radius. Increasing Rcfe modestly increases average source radius at all times, and the value of A only affects source radius at later 
times when the disk is depleted of chondrules. 


sity (due to smaller diffusion length) can be envisioned. 
This would be compatible with a temporal evolution as well 
( [Jacquet, Gounelle &: Fromang|2012p . The observed preser¬ 
vation of refractory inclusions, which formed early in the so¬ 


lar system (perhaps within the first 0.1 Myr; e.g. Bizzarro, 
[Baker &; Haack||2004[ [Amelin et al.|[2Q10p , in large abun¬ 
dances in CCs may be a further argument in favour of low a 
values i Jacquet, Fromang fc Gounelle|2011 K Still, it should 
be noted that simulations of refractory inclusion transport 
by Yang & Ciesla (20121 achieve suitable levels of preser¬ 
vation assuming a moderately high a = 10“^, provided the 
disk was initially very compact {R\ < 10 AU) so that its vis¬ 
cous expansion sent many refractory inclusions far from the 
Sun (10-100 AU), hence a longer drift timescale afterward. 
In the current understanding of disk magnetohydrodynam¬ 
ics, a dead zone would nonetheless have rapidly emerged 
anyway, strengthening the reported retention. 

A further noteworthy constraint is provided by analyses 
of dust returned from comet Wild 2 by the Stardust mission. 
Indeed, while a carbonaceous chondrite-like component is 
certainly present, oxygen isotopes (Nakashima et al. 20121 
and chemical compositions ( [Frank, Zolensky &; Le||2014p of 
olivine grains show that the contribution of most carbona¬ 
ceous chondrite groups is limited, and that CR chondrites 
as well as non-carbonaceous groups (the former showing 
some characteristics transitional with the latter; j Jacquet 


Robert 20131 may be significant sources. This is at vari¬ 


ance with the idea that non-carbonaceous and carbonaceous 
chondrites only differ by their spatial origins, with carbona¬ 
ceous chondrites being furthest, since the spatial source dis¬ 
tribution should converge in the outer solar system (see Sec¬ 
tion 4.21, so that in terms of CFR-originated components. 


comets should be dominated by carbonaceous chondrite 
ones. The Stardust mission may thus offer further evidence 
for a temporal evolution, with possibly Wild 2 being the re¬ 
sult of a relatively late accretion (Nakashima et al. 20121, 
if the suggestion by Jacquet, Gounelle & Fromang (20121 
that noncarbonaceous chondrites accreted later than their 
carbonaceous counterparts holds. This may also explain the 
lack of evidence of initial live ^®A1 for several Wild 2 parti¬ 
cles (e.g. Ogliore et al.|20l2 l. Carbonaceous chondritic ma¬ 
terial may have survived longer at greater distance, e.g. in 


the wake of an initial expansion of the disk (Yang & Ciesla 


20121 . 


6 CONCLUSIONS 

We have performed a numerical investigation of chondrule 
transport in the young Solar System. Our results speak to 
the nature of the chondrule-matrix relationship; spatial and 
temporal constraints on the chondrule formation process; 
and the success of deploying astrophysical numerical meth¬ 
ods for cosmochemical applications. The main results are as 
follows: 

(i) We have defined a “complementarity parameter” (^) 
as a metric for comparing the chondrule-matrix relationships 
in our astrophysical models to those found in lab-analysed 
chondrites. It evaluates how close the chondrule and ma¬ 
trix component are to being complementary, while being 
independent of how close the individual components are to 
having solar abundances. 

(ii) Our simulations showed that the gas-solid decoupling 
parameter S (the ratio of the diffusion to drag timescales) is 
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predominantly responsible for these relationships, with suf¬ 
ficiently low values of S (the diffusion dominated regime) al¬ 
lowing for the time between chondrule production and chon- 
drite accretion to be longer than the disk’s viscous timescale. 
Older disks, with lower mass accretion rates, are heavily de¬ 
pleted of chondrules relative to their younger counterparts. 
Thus, we constrain the chondrule formation events to have 
occurred early, in a disk with M > 10“® Mq yr“^ for nom¬ 
inal parameters. 

(hi) At a given radius and time in the planetary disk, 
chondrule origin varies as a function of space and time, each 
of which are affected by disk parameters as well as the rate 
of chondrule formation. Again, mass accretion rate, and thus 
disk age, is seen to be relevant to chondrule distributions. 

(iv) The location of the S = 1 line limits the outward 
diffusion of chondrule matter. The distribution and extent 
of chondrule formation in space significantly impacts spa- 
tiotemporal diversity, because most chondrules accreted in 
the inner disk formed locally. 

(v) The spatial diversity in our simulations is insufficient 
to explain the number of distinct chondrite classes that ex¬ 
ist, but we argue that temporal diversity can, and likely 
does, account for this discrepancy. More radiochronological 
data establishing temporal bounds on chondrite classes in 
terms of either chondrule age or chondrite accretion time 
are obviously needed. 

This study shows the viability of our numerical meth¬ 
ods in investigating chondrule origins. The agreement be¬ 
tween our simulations for a single CFE, continuous CFEs 
with a static disk, and continuous CFEs with an evolving 
disk convey the power of this technique for exploring regimes 
unattainable through pure analytics or cosmochemical spec¬ 
ulation, and we note that such simulations are not computa¬ 
tionally expensive. Furthermore, our simulations make pre¬ 
dictions of their own, which can be candidly evaluated in 
future studies. This successful marriage of a cosmochemical 
quandary to an astrophysical method shows the strength of 
this relationship, and we ardently promote the exploitation 
of similar relationships in future endeavours. 
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APPENDIX A: CALCULATION OF Thalf IN 
THE REGIME 1 < S' < Rcentre/L 

In this regime (as can be verihed a posteriori), the initial 
widths of the chondrule and dust populations can be ne¬ 
glected and the variation of heliocentric distance can be ig¬ 
nored. After a time t, the surface density profiles of chon- 
drules and dust can be approximated by a Gaussian distri¬ 
bution: 

rr>\ Mch,d ( “ -Rcentre,-i,,d) /.in 

= U--j- 


where a = V2Dt for both populations, and Mch and Md are 
the masses of the chondrule and dust populations, respec¬ 
tively. The complementarity parameter then reduces to: 

>. ^ch/^d / “h d (RcentrefU -^centrej) 

- 

(A2) 

with 


d — Rcentred Rcentre,,!, — (Udriftd ^drifted ) ^ 


~ I'Udriftei, 1^- 


(A3) 


from which one can obtain the abscissas corresponding to 
^ = a and = b. Integrating then the chondrule surface 
density in between yields 


orf I £. — 1 .4 1 _ prf ( llLS; 2 _ 1 d \ 

^ \V2 d 2vT n/ 2 d 2 V 2 J 

WaS^C^t, — -X-, 

(A4) 

giving I « 2.297 for {w,a,b) = (0.5, 0.7,1.3). Plugging thaif 
into the definitions of a and d and substituting for Udrift^h 
gives 


2D ( P 
2.2972 \T-dP/dR 


2 


(A5) 
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